On the connection between structural distortion and magnetism in graphene with a 

single vacancy 
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The correlation between structural distortion and emergence of magnetism in graphene containing 
a single vacancy was investigated using first-principles calculations based on density functional 
theory (DFT). Our results have shown that a local distortion is formed around the vacancy, with 
reconstruction of two atomic bonds and with a dangling bond remaining at the third atom adjacent 
to the vacancy. A systematic investigation of the possible out-of-plane displacement of this third 
atom was then carried out, in order to ascertain its effects on the magnetic features of the system. 
The ground state was definitely found to be magnetic and planar, with spin-resolved a and n bands 
contributing to the total magnetic moment. However, we have also found that metastable solutions 
can be achieved if an initial shift of the third atom above a minimum threshold from the graphene 
plane is provided, which leads to a non-planar geometry and a non-magnetic state. 



T3 

o 
o 



>: 

in, 
ov 



(N 



X 



I. INTRODUCTION 



The issue of magnetism in graphene and re- 
lated materials-including graphene multilayers, 
graphite, nanocarbons and others-has been the fo- 
cus of intense research in the last years, both from the 
experimental and the theoretical points of view [H-Q- 
The onset of intrinsic magnetic properties is always 
linked to some kind of defects introduced in the bipartite 
lattice, such as atomic vacancies, chemisorbed species or 
edge sites 0, H, @ . Single vacancies cause the appearance 
of dangling bonds and, depending on the defect con- 
centration and on the degree of passivation associated 
with eventually present hydrogen atoms, the magnetic 
moments due to electrons with uncompensated spins in 
a and tt orbitals can couple to give a total magnetic 
moment in the range 1-2 /ib per atomic vacancy [H,H, 0] ■ 
Although there have been some experimental reports 
suggesting the existence of ferromagnetism at room 
temperature in defective graphene [a Q , recent work 
using measurements of magnetic properties [Toj and 
muon spin rotation spectroscopy has ruled out this 
possibility. In particular, Nair et al. showed that 
localized defects carry magnetic moments with spin i, 
leading to strong paramagnetism but no magnetic order 
down to liquid He temperature. Some recent calculations 
performed by Palacios & Yndurain for variable defect 
concentrations in graphene sheets containing single 
vacancies have confirmed these findings, by showing 
that the vacancy-induced magnetism associated with 
extended tt orbitals has a trend to vanish at any exper- 
imentally relevant vacancy concentration-thus spoiling 
any possibility of ferromagnetic or antiferromagnetic 
interactions associated with those tt electrons Q. On 
the other hand, this work also showed that localized, 
non-interacting states associated with a flat a band 
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are the only ones surviving for defect concentrations 
approaching zero, in apparent agreement with the 
paramagnetic properties experimentally observed for 
graphene irradiated with protons [l(| • 

When a single atomic vacancy is formed in an oth- 
erwise perfect two-dimensional graphene lattice in the 
absence of hydrogen (or other gases that could lead to 
passivation), three dangling bonds remain in the atoms 
surrounding the vacancy. A Jahn- Teller distortion in this 
neighbourhood leads to the reconstruction of the bond 
between two of the atoms closest to the vacancy, leav- 
ing the third atom less bound to the network and with 
an unsaturated dangling bond [H, [HI], as illustrated in 
Figure Q] The question then arises whether or not this 
atom (identified as atom 3 in Figure 1) will remain in 
the plane, an issue that has received a number of appar- 
ently conflicting responses in the last years. El-Barbary 
et al. used the local density approximation (LDA) in an 
atomic cluster to model a defective hydrogen-terminated 
graphene ribbon [l3| . They found that the fully op- 
timized structure is distorted, with the atom with the 
unsaturated dangling bond moving out of the plane by 
0.47 A. When using calculations with spin polarization, 
the spin-polarized state was found to be higher in en- 
ergy than the unpolarized one by ca. 0.5 eV, leading the 
authors to conclude that the ground state of vacancy- 
containing graphite should be nonmagnetic (T3 |. A simi- 
lar method was followed by Dharma-wardana & Zgierski, 
who used DFT calculations to study the structural, elec- 
tronic and magnetic properties of a single vacancy in a 
finite graphene fragment with zig-zag edges terminated 
with hydrogen atoms |14(. They found the ground state 
was spin polarized and planar, lying ca. 0.2 eV below the 
non-planar unpolarized structure, with the latter show- 
ing a saddle-like distortion around the vacancy. In con- 
trast, when using periodic boundary conditions to study 
an infinite graphene sheet containing a single vacancy, 
no out-of-plane distortions were found in the fully re- 
laxed structure, no matter whether or not spin polariza- 
tion was used in the calculations [l4| . In a related work, 
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Ma et al., employing the generalized gradient approxima- 
tion (GGA), concluded that the ground state of graphene 
containing a single vacancy was spin polarized, with the 
atom with no reconstructed bonding displaced by 0.18 
Aout of the plane [l5j . When forcing a nonmagnetic so- 
lution, which was found to be higher in energy by 0.1 
eV as compared to the ground state, they found this dis- 
placement being increased to 0.46 A, in agreement with 
the results of El-Barbary et al. [l3j . Accordingly, Dai et 
al. also found, using the same GGA method, a magnetic 
moment of 1.33 /i B for a monovacancy in graphene and 
reported an upward shift of atom 3 from the plane by 
0.184 A[ll|. These results are in contrast with the very 
recent GGA calculations reported by Nanda et al., who 
found no displacement of any atom out of the graphene 
plane containing a single vacancy [l7j . 

A more detailed investigation of this issue was reported 
by Faccio et al., who used GGA calculations to com- 
pare the effects caused by carbon vacancies and boron 
doping in graphene [18] . By studying monovacancies in 
graphene sheets of different sizes (and thus with variable 
degree of interaction between vacancies in different super- 
cells), these authors found that a nonmagnetic solution 
was higher in energy by ca. 0.4 eV as compared to the 
magnetic structure for the whole range of defect concen- 
trations. Also, the establishment of the nonmagnetic so- 
lution was accompanied by an out-of-plane displacement 
of atom 3 by 0.3 A. Interestingly, it was reported that 
when this nonmagnetic solution was used as the starting 
point of a new spin-polarized calculation, allowing full 
structural relaxation, then the structure relaxed back to 
the flat graphene sheet, with a net magnetic moment in 
the range 1.0-1.3 (ib, depending on the defect concentra- 
tion [l8[. The results recently reported by Ugeda et al. 
in an investigation about the effects of single vacancies in 
a graphene monolayer deposited on Pt surface [lj| also 
helped to shed some light into this subject. In order to 
understand the absence of magnetism in the graphene/Pt 
system, the authors first studied the features of a pure 
graphene sheet containing a single vacancy; again, spin- 
polarized calculations showed that the ground state was 
magnetic (with a magnetic moment of 1.5 /is) and nearly 
planar, with a nonmagnetic solution found 0.1 eV above 
the ground state. However, by systematically forcing an 
out-of-plane displacement of atom 3, a decrease in the 
magnetic moment as a function of the displacement was 
obtained, with a nonmagnetic state being found (even 
using spin-polarized calculations) for a displacement of 
ca. 0.5 A. This quenching of magnetic moment was at- 
tributed to the mixing of a and 7r states in the spin- 
resolved density of states (DOS), causing a partial or 
total compensation of the magnetic moments associated 
with the former a and 7r orbitals [19|. This conclusion 
was in agreement with the work by Dharma-wardana & 
Zgierski, who reported that in the unpolarized structure 
of a finite graphene fragment with a single vacancy the 
distortion caused a disconnection of the electron associ- 
ated with the vacancy from the remaining 2D electron 



network; this led to the formation of two singlet electron 
pairs - one associated with the two electrons at atom 
3 and the other one associated with the reconstructed 
bonding between atoms 1 and 2 (see Fig. [1]), causing the 
structure to become nonmagnetic [14j . 

Considering this scenario with some apparently dis- 
crepant conclusions as well as a number of sparse comple- 
mentary findings, we decided to undertake a systematic 
investigation of the effects of out-of-plane displacements 
on the structural, electronic and magnetic features of a 
graphene sheet containing a single vacancy. Both spin- 
polarized and unpolarized calculations with full struc- 
tural relaxation were performed, allowing us to compare 
the total energies, magnetic moments and electronic con- 
figurations of each relaxed structure corresponding to dif- 
ferent initial displacements. These results unequivocally 
show that the ground state is indeed planar and mag- 
netic, but the presence of other local energy minima can 
lead to metastable solutions with different structural and 
electronic properties. 



II. COMPUTATIONAL DETAILS 

In this work, we have investigated carbon vacancies 
in a graphene sheet, through first-principles^ calculations, 
based on density functional theory (DFT)[1|. The DFT 
calculations were performed using ultrasoft Vanderbilt 
pseudopotentials [2l| , and a generalized gradient approx- 
imation (GGA) for the exchange-correlation potential 
[H, as implemented in the VASP code dlH. In order 
to study the vacancies we have used 6x6 supercell with 
72 carbon atoms. The vacancies were created by remov- 
ing one carbon atom in the graphene sheet. The lattice 
parameter for graphene obtained from structural opti- 
mization was 2.46 A. We have used a plane-wave-cutoff 
energy of 400 eV and a Monkhorst-Pack [26[ scheme with 
a 5x5x1 k-mesh for the Brillouin zone integration. Cal- 
culations were performed both with and without spin po- 
larization for the graphene sheet containing a single va- 
cancy. In all calculations the lattice parameter was kept 
fixed at the calculated value, whereas the atoms were al- 
lowed to relax until the atomic forces were smaller than 
0.025 eV/A. 

III. RESULTS AND DISCUSSIONS 

Fig[TJ shows the fully relaxed atomic structure of the 
graphene sheet with a single carbon vacancy, obtained 
by spin-resolved DFTcalculations. 

This structure clearly shows the local breaking down of 
the threefold symmetry due to the Jahn- Teller distortion 
induced by the reconstruction of two dangling bonds left 
after removing one carbon atom from the graphene sheet. 
It should be noted that, as a result of this distortion, the 
atoms 1 and 2 move closer to one another, forming a 
reconstructed bond with length of 2.07 A(as compared 
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FIG. 1: The relaxed atomic structure of the graphene sheet 
with a single vacancy, where spheres represent the carbons 
atoms. The arrow indicates the atom 3 with the dangling a 
bond and atoms 1 and 2 are those rebonded. 



to 2.46 A in defect-free graphene). This gives rise to 
the in-plane displacements of other carbon atoms in the 
relaxed structure, as shown in Table HI 



TABLE I: Distances corresponding to first and second neigh- 
bors in a graphene sheet with a single vacancy and in defect- 
free graphene, where di-j denotes the distances corresponding 
to different pairs of atoms as shown in Figffl 



Graphene-vacancy (A) Graphene (A) 



First-neighbors 



d 3 _ 4 =d3-i2 =1.37 
d 2 -7=di_ 9 =1.40 
d 8 -7=d8-9 =1.41 
d5_4=dll-12=da-6 = 1.42 
d6-5=dl0-ll=da-c = l-43 
d2-6=rfl-9=rfl0-l =1.47 



do = 1.42 



Second- neighbors 



di_ 2 =2.07 
d 7 -9=2.37 
d2-8=e£8-i=2.39 
d 5 -3=dn-3=2.43 
d 6 -4=dio-i2=2.44 
di-3=e£2-3=2.56 
d 2 -s=2.58 
dn_i=2.59 



d =2.46 



Table U shows that there exists a variation in the in- 
teratomic distances corresponding to second neighbors, 
ranging from 2.07 A up to 2.59 A . The largest change in 
the in-plane interatomic distances is observed for atoms 
1 and 2 (di_2=2.07 A) and is associated with the recon- 
structed bond mentioned above. We can also note that 
the interatomic distance between first neighbors varies 
from 1.37 A up to 1.47 A , compared to 1.42 A in defect- 
free graphene. These results are in good agreement with 
previous work 0, [HI ■ 



TABLE II: Final (Z'/ and Z/ ) distances of atom 3 measured 
from the sheet after full relaxation for unpolarized and spin- 
polarized calculations, respectively. (Zi) represents the initial 
displacement of atom 3 from the graphene sheet and (/i) is 
the magnetic moment corresponding to the relaxed structure. 

Z, (A) Wf % (A) \ (pL B ) 



0.00 


0.00 


0.00 


1.21 


0.06 


0.40 


0.02 


1.21 


0.08 


0.40 


0.05 


1.20 


0.10 


0.41 


0.06 


1.19 


0.20 


0.41 


0.11 


1.16 


0.30 


0.41 


0.13 


1.12 


0.40 


0.41 


0.20 


1.04 


0.50 


0.42 


0.41 


0.00 


0.60 


0.42 


0.41 


0.00 


0.70 


0.43 


0.42 


0.00 



The formation energy of a single vacancy in the 
graphene sheet at the ground state of the system was 
calculated as 



E 



f 



-{E v 



N -n 
N 



Eg), 



(1) 



where N denotes the number of C atoms in the defect- 
free graphene sheet, n represents the number of vacant 
C atoms (i.e., n = 1 in the present case) and E g and E v 
are the total energies of the defect-free and the vacancy- 
containing graphene sheets, respectively. The formation 
energy of the graphene sheet containing a single vacancy 
calculated according to eq.fll]) was found to be 7.6 eV, 
which is in good agreement with the experimental value 
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of 7.0 eV 27 1 and also with previous results of DFT cal- 



culations 

In order to verify the existence of possible metastable 
solutions involving different geometries of the graphene 
sheet containing a single vacancy, we have investigated 
in detail the possible displacement of atom 3 (which is 
the atom that does not form a reconstructed bond with 
the other ones, indicated by the arrow in Fig. Q]) per- 
pendicularly to the sheet. Then, a series of calculations 
were performed, starting from a structure with atom 3 
initially displaced from the sheet by variable distances 
(named Zj). These calculations were done for both un- 
polarized and spin-polarized schemes, with the structure 
being allowed to fully relax. The initial (Zj) and final po- 
sitions (named Z'/ and Z/ for the unpolarized and spin- 
polarized calculations, respectively) of atom 3 with the 
corresponding magnetic moments are summarized in Ta- 
ble HH whereas the plots of the total energy of the system 
as a function of the final position reached by atom 3 are 
shown in Fig. [2] 

From the data shown in Table [TT| we can analyze the 
behavior of the final position reached by atom 3 as a 
function of its initial position. In the unpolarized case, we 
can observe that if an initial out-of-plane shift is given to 
atom 3, then its final position reaches values around 0.42 
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A(even for initial values as small as 0.06 A). On the other 
hand, for spin-resolved calculations there exist different 
metastable solutions that can be reached depending on 
the initial shift. In particular, it should be noted that 
for Zj < 0.40 A the final position of atom 3 Zj varies 
from 0.0 up to 0.20 A , whereas for Zj > 0.40 the final 
position is around 0.42 A, matching the value found in 
the unpolarized calculations . 

Fig. [5] shows how the total energy of each relaxed struc- 
ture changes as a function of Z'/ or Zf, for calculations 
performed without (Figure [2^) or with (Figure^) spin 
polarization, respectively. These results show that the 
ground state of the relaxed structure is indeed planar 
(i.e., Zf = 0) and spin-polarized, with a total energy 
nearly 0.10 eV below that of the unpolarized structure, 
in good agreement with some previous results [lEI . [l^ | . 
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the ground state), which can help to understand the re- 
ports of shifts with this magnitude found in the literature 
[T3L [l5| . Further, it can be observed in Figure [5b that 
a little energy variation and small differences between 
Zi and Zf are detected among the several final struc- 
tures obtained for Zj in the range from 0.00 up to 0.20 
A , which is also consistent with other reports [HI, [l6[ . 
These results show that the final reached structures are 
strongly dependent on the use of spin polarization in the 
calculations [H, and that, even when using spin po- 
larization, the setting of the initial conditions can signif- 
icantly change the properties of the relaxed structure. It 
is likely that most of the discrepancies found in the liter- 
ature regarding the out-of-plane displacements of atom 3 
can be associated with the trends illustrated in Fig. [3J 

Moreover, the abrupt change of the energy for Zf > 
0.20 A can be associated with the reported disconnec- 
tion of the 7r electron associated with the vacancy from 
the 2D electron network [l4j]. For such large shifts, it 
is energetically more favorable for the system to accom- 
modate two of the electrons released by the vacancy in 
a singlet pair at atom 3, with the other two electrons 
becoming part of the reconstructed cr-like bond between 
atoms 1 and 2. Thus, this metastable structure must 
remain nonmagnetic, with the same energy (-650.16 eV) 
irrespective of the use of spin polarization in the calcula- 
tions. This behavior can be clearly observed in Figure 02 
where the magnetic moment of the graphene sheet con- 
taining a single vacancy is plotted as a function of Z/ . 
It is visible that the planar ground state (Z/ = 0.00 A) 
is magnetic, with a magnetic moment of 1.21 /i^; on the 
other hand, for Zf > 0.20 A the system evolves to a 
nonmagnetic state, achieved with Zf close to 0.42 A. 



FIG. 2: Total energy of the graphene sheet containing a single 
vacancy calculated as a function of the displacement of atom 
3 perpendicular to the sheet (a) without spin polarization and 
(b) with spin polarization. The lines connecting between the 
points are guides for the eyes. 

As was already mentioned before, the details of the re- 
laxed structure of a graphene sheet containing a single va- 
cancy reported in the literature vary widely. While some 
authors have reported a planar geometry for the ground 
state there are also many reports of non-planar 

structures with out-of- plan e displacements varying from 
0.18 to 0.47 A [HElll!]. lt is worth emphasizing that 
in the case of the unpolarized calculations (Figure [3^) , 
atom 3 tends to be shifted by 0.42 A out-of-plane, 
whereas the spin-polarized calculations (Figure [5b) yield 
a planar structure. However, it is interesting to note that 
even in the case of spin-polarized calculations there exists 
a clustering of points close to this value of Zf, suggesting 
a local energy minimum in this region. This means that 
if the initial shift of atom 3 is above ca. 0.20 A , the 
structural relaxation leads to metastable states with Z f 
around 0.42 A and with the same total energy as com- 
pared to the unpolarized structure (ca. 0.10 eV above 
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FIG. 3: Magnetic moment of the graphene sheet containing a 
single vacancy calculated as a function of the displacement of 
atom 3 perpendicular to the sheet for the relaxed structure. 
The line connecting the points is a guide for the eyes. 

The changes in the electronic structure related to dif- 
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ferent atomic configurations as atom 3 is moved away 
from the graphene sheet are shown in Fig. 21 where we 
can see the spin resolved band structure for the different 
Zf values, after full relaxation of the system. 




r m k rrMK r r m k r 



FIG. 4: Spin-resolved band structures calculated after full 
relaxation of the graphene sheet containing a single vacancy 
for three different values of the final displacement of atom 3 
perpendicular to the sheet: (a) Zf = 0.00 A ; (b)Z/ = 0.20 
A and (c) Z f = 0.42 A. The zero of energy is at the Fermi 
level. Solid and dotted lines represent spin-majority and spin- 
minority bands, respectively. 

The existence of separated spin-minority and spin- 
majority bands is evident in Fig. 0Ji and, to a lesser 
extent, in Fig. |4Jd, both for a and n bands, corresponding 
to Z f = 0.0 A and 0.20 A, respectively. With increasing 
Zf, one can also observe a broadening of the a band that 
was located close to -0.5 eV for Zf =0 and that was as- 
sociated with a strongly localized state around atom 3 in 
the planar structure []], Gil ■ At the same time, the spin- 
minority 7r band is pushed down and both 7r bands are 
almost completely below the Fermi level for Zf = 0.20 
A, which leads to a reduction in the magnetic moment 
compared to the planar geometry, as shown in Figure [3] 
Moreover, Fig. 3J: shows that for Zf — 0.42 Athere exists 
a complete overlap of spin-majority and spin-minority 
bands, leading to the quenching of the magnetic moment 
of the sheet. Then, one can note the a and ir states ob- 
tained for Zf <0.20A evolved to a* and ir* states, which 
are combinations of bonding and antibonding orbitals. 
This leads to a stable electronic configuration with dou- 
ble occupancy of the a* state and, partially, of the tt* 
states, which are now pinned to the Fermi level. The 
quenching of the magnetic moment of the sheet associ- 
ated with the out-of-plane displacement of atom 3 has 
already been previously reported [lj|. However, it is not 
clear whether or not a full relaxation of both atomic and 



electronic coordinates was allowed in the structures re- 
ported in ref. [19( , so it is difficult to quantitatively com- 
pare those results with the ones presented here. But the 
qualitative mechanism describing the quenching of the 
magnetic moment is clearly consistent with the trends 
shown in Figs. OE]and[I] as atom 3 is shifted away from 
the graphene plane, there exists a progressive change of 
hybridization of that atom and a mixture of a and tt 
states; this changes the occupation of the corresponding 
spin-majority and spin-minority bands and thus leads to 
a partial compensation of the spin populations and a re- 
duction in the magnetic moment of the sheet. In the 
present case, when the displacement of atom 3 exceeds a 
threshold close to Zf = 0.20 A , then a complete com- 
pensation of spin populations in the now mixed states 
occurs, leading to the settling of a nonmagnetic state. 



IV. CONCLUSION 

DFT calculations have been employed in a systematic 
investigation of the effects of out-of-plane atomic dis- 
placements on the structural, electronic and magnetic 
features of a graphene sheet containing a single vacancy. 
The results showed the occurrence of a local distortion 
around the vacancy in the relaxed structure, with recon- 
structed atomic bonds between two atoms close to the 
vacancy and with the third atom located in the graphene 
plane. Total energies, magnetic moments and electronic 
configurations of the final structures reached after full re- 
laxation starting from different initial shifts of this third 
atom in the direction perpendicular to the graphene sheet 
were calculated, both for polarized and also for unpo- 
larized structures. The present results showed that the 
ground state of this system is indeed planar and mag- 
netic, with spin-resolved a and tt bands contributing to 
the total magnetic moment of the graphene sheet con- 
taining a single vacancy. However, we have found the 
presence of other local energy minima, which can lead to 
metastable solutions with different structural, electronic 
and magnetic properties that are strongly dependent on 
the extent of the out-of-plane shift. 
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